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We study a dilute mixture of degenerate bosons and fermions across a Feshbach resonance of 
the Fermi- Fermi scattering length ay. This scattering length is renormalized by the boson-induced 
interaction between fermions and its value is crucial to determine the phase diagram of the system. 
For the mixture in a box and a positive Bose-Fermi scattering length, we show that there are three 
possibilities: a single uniform mixed phase, a purely fermionic phase coexisting with a mixed phase, 
and a purely fermionic phase coexisting with a purely bosonic one. As 1/clf is increased from a 
negative value to the Feshbach resonance (1/ap = 0) the region of pure separation increases and 
the other two regions are strongly reduced. Above the Feshbach resonance (1/ap > 0), pairs of 
Fermi atoms become Bose-condensed molecules. We find that these molecules are fully spatially 
separated from the bosonic atoms when 1/of exceedes a critical value. For a negative Bose-Fermi 
scattering length we deduce the condition for collapse, which coincides with the onset of dynamical 
instability of the fully mixed phase. We consider also the mixture in a harmonic trap and determine 
the conditions for partial demixing, full demixing and collapse. The experimental implications of 
our results are investigated by analyzing mixtures of 6 Li- 23 Na and 40 K- 87 Rb atoms. 

PACS numbers: PACS Numbers: 03.75.Ss, 03.75.Hh, 64.75. +g 



I. INTRODUCTION 

The regime of deep Fermi degeneracy is now actively 
studied with ultracold vapors of 6 Li and 40 K atoms. Ex- 
periments on two-hyperfine-state Fermi gases [l|, H H H| 
are concentrated across a Feshbach resonance, where a 
crossover from a Bardeen-Cooper-Schrieffer (BCS) super- 
fluid to a Bose-Einstein condensate (BEC) of molecular 
pairs has been predicted H H @] . Experimental and the- 
oretical investigations of the Fermi cloud across the BCS- 
BEC crossover have been devoted to density profiles [1 H 
OH , collective excitations M IIEO, GJ, El Bill EpLcon- 
densate fraction fi^, El fl9ll20||. free expansion UH], Hl[ 
and vortices [22|, [23| . 

An interesting issue is the inclusion of Bose atoms, 
like 23 Na or 87 Rb isotopes, in the Fermi cloud. Trapped 
boson-fermion mixtures, with Fermi atoms in a single hy- 
pcrfine state, have been investigated by various authors 
both theoretically M, [p, [H M, HI, M, H M, HI and 
experimentally [33|, [34|, [35| . Very recently McNamara et 
al. [3t| have reported the observation of simultaneous 
quantum degeneracy in a completely spin polarized di- 
lute gaseous mixture of 3 He and 4 He in their first excited 
metastable states [37| 

The purpose of this paper is to carry out a study of 
the miscibility of a Fermi-Bose mixture with the Fermi 
atoms in two equally populated hyperfine states across a 
Feshbach resonance of the Fermi-Fermi scattering length. 
Since the density of spin up and spin down fermions is 
bound to be the same, this three-component system be- 
haves effectively as a two-component mixture: a super- 
fluid Fermi gas and a Bose-Einstein condensate. 

First we consider the mixture in a box. For a positive 
Bose-Fermi scattering length we show that, depending 
on the total number densities rip and ns of the Fermi 



and Bose components, there are three possible minimum 
energy configurations: a uniform phase with bosons and 
fermions fully mixed, a purely fermionic phase coexisting 
with a purely bosonic one, and a purely fermionic phase 
coexisting with a mixed phase. We find that the region in 
the hf - tlb plane where each one of the three is stable, 
strongly depends on the effective Fermi-Fermi scattering 
length. In addition, we find that when the Bose-Fermi 
scattering length is negative, collapse is driven by the 
dynamical instability of the homogeneous mixture. 

We discuss also the mixture under harmonic confine- 
ment. In this case the conditions for partial demixing, 
full demixing and collapse crucially depend on the sign 
of the Bose-Fermi scattering length, the ratio between the 
Bose-Fermi and the Bose-Bose scattering lengths, and the 
number of atoms involved. We analyze the experimen- 
tal implications of our results by considering mixtures of 
6 Li- 23 Na and 40 K- 87 Rb atoms as examples. 



II. RENORMALIZED SCATTERING LENGTH 

It has been recently shown [lj| El [H, H^] that at zero 
temperature the energy density of a uniform and dilute 
Fermi superfluid, composed by two equally populated 
spin states, can be written in the BCS-BEC crossover 
as 
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AnTf{y) , 



(1) 



where A = h 2 (3-7T 2 ) 2 / 3 / (2m f) and mp is the mass of a 
fermionic atom. The universal function f(y) depends on 
the inverse interaction parameter y — (kp op) , where 
ap is the Fermi- Fermi scattering length, uf is the Fermi 
number density and kp = (S^np) 1 / 3 is the Fermi wave 
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vector. In the experiments on the BCS-BEC crossover 
the scattering length a F is changed by using an exter- 
nal magnetic field (Feshbach resonance technique) and 
may be varied from large negative to large positive val- 
ues 0, 0, @, 0] ■ In these experiments the effective range 
of the interaction between fermions is much smaller than 
the mean interparticle distance and so the Fermi cloud 
may be considered as dilute. In the weakly attractive 
regime (y <C —1) there is a BCS Fermi gas of weakly 
bound Cooper pairs and the universal function f{y) has 
the asymptotic behavior f(y) = 1 + 10/(97ry) + 0{l/y 2 ) 
[3 EE Hl| ■ In the so-called unitarity limit (y = 0) one 
expects that the energy per particle is proportional to 
that of a non-interacting Fermi gas with a coefficient 
/(0) = 0.42 [1(3]. In the weak-coupling BEC regime 
(y ^> 1), there is instead a weakly repulsive Bose gas of 
molecules of mass tum — 2m F and density um — n F /2. 
Such Bose-condensed molecules interact with a positive 
scattering length ajvf = 0.6a F [38l . l4l| and the univer- 
sal function f(y) has the asympotic behavior f(y) — 
5a M /(18na F y) + 0(l/y 5 / 2 ) (jl, OJ, [H|. Actually the 
formula bm = 0.6a F has been analytically derived for a 
weakly-interacting gas of dimers [4lj but it seems to work 
well also in the strongly-interacting limit [38| . 

In the present work we investigate the effect of Bose 
atoms mixed with the two-hyperfine-state gas of fermions 
by using the two different analytical forms Eq. ^ and 
Eq. ([3]) of the universal function f(y) proposed by Kim 
and Zubarev [l5| and by Manini and Salasnich [lgj to fit 
the numerical results from Monte Carlo (MC) simulations 
[H, Hi| and obeying the above asymptotic expressions. 

The formula proposed by Kim and Zubarev is a [2/2] 
Pade approximant 



f(y) = do 



diy + d 2 



u 



d^y + d 4 



(2) 



where the values of the five parameters di are reported 
in Ref. [l5T| , Manini and Salasnich slightly improved the 
fit to the available MC data [H, Hi| by using a different 
parametrization, i.e. 



f(y) = ai - a 2 arctan a 3 y 
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with the values of the parameters and /3j reported in 
Ref. 0. 

At zero temperature a dilute gas of bosons, with den- 
sity n B , mass mj and interaction strength g B , is a Bosc- 
Einstein condensate, with a coherence length £b given 
by £b = fr/(2m B c B ), where cb — {g B n B /tub) 1 ^ 2 is the 
bosonic sound velocity. Note that the strength g B is given 
by gB = 4:Trh 2 a B /m B , where as is the Bose-Bose scat- 
tering length. Already back in the 1960's, in connection 
with dilute solutions of 3 He in superfluid 4 He, it was real- 
ized that [13] the Fermi-Fermi scattering length is renor- 
malized by the density fluctuations of the bosons and that 
this renormalization may lead to a BCS state even when 
the Fermi-Fermi bare interaction is repulsive. A thor- 
ough analysis [HI, [27, 28, 29] of the renormalization of the 
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FIG. 1: Inverse renormalized scattering length l/a' F as a func- 
tion of the inverse bare scattering length l/a F , according to 
Eq. (UJ). In the inset there is linear plot of a' F as a function 
of of- Lengths are in units of g% F / '(Ags)- 



Fermi-Fermi scattering length in Fermi-Bose mixtures of 
ultracold atomic gases has been provided in recent years. 
Viverit 29] has shown that, quite remarkably, in the di- 
lute mixture and under the condition k F £,B 1, the 
renormalized scattering length of the Fermi-Fermi inter- 
atomic potential does not depend on the densities of the 
atomic species but only on their interaction strengths ac- 
cording to 



a F = a F 



m F g BF 
Anh 2 gB 



(4) 



where gsF = 47r?i asF I^bf is the Bose- Fermi inter- 
action strength, ci B f its scattering length and itibf = 
2m B mF/(m B +m F ) the reduced mass. Heaviside's func- 
tion Q(rib) merely states that the renormalization occurs 
only where / 0. 

In Fig. 1 we plot the renormalized l/a' F as a func- 
tion of l/a F according to Eq. Fig. 1 shows that 
the Bose-Fermi interaction strength gsF drives the sys- 
tem into the non- interacting regime (l/a' F — ±oo) when 
= 9 2 B F^F/(9B^h 2 ) = {%/{^y/3){Ag B /g 2 BF ). 
As previously stressed, the renormalization formula 
holds if k F ^B <C 1, which may be satisfied only when g B 
is positive since it is equivalent to 



n F <€. 



f 2m B g B \ 
\ m F A J 
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(5) 



In the following we will use adimensional densities: for 
fermions n F — n F g BF /(A 3 g B ) and for bosons n B = 

n B g BF /(A 3 g B ), so that Eq. ([5|) becomes n F <C an 3 B 2 , 
where a = (2m B g B /(m F \g BF \)) 3 ^ 2 . In the case of a mix- 
ture of 6 Li and 23 Na atoms, where g BF /g B — 0.84 [35| . 
one finds a ~ 27, while for a mixture of 40 K and 87 Rb 
atoms, where g BF /g B — —4.56 [H[, one gets a ~ 1. In 
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both cases the range of densities where this inequality is 
valid is sufficiently large to allow experimental tests on 
the predictions of the following sections. Note that mea- 
suring lengths in units of g%p/(AgB), Eq. (4) becomes 
a F = ap — (97r) 1 / 3 /8 ~ ap — 0.38; working near the Fes- 
hbach resonance (a' F — cif — oo) the renormalization of 
the Fermi-Fermi scattering length is very small. 

The Eq. (0| is valid for both negative and positive 
ap under the condition ( [5t pr ovided that a Fermi- Fermi 
bound state is not formed [2!| • We expect that the effect 
of renormalization is negligible close to the resonance also 
on the BEC side of the BCS-BEC crossover, where indeed 
there is the molecular bound state made of two Fermi 
atoms in different hyperfine states. Thus, studying the 
BCS-BEC crossover, in the BEC side (ap > 0) we take 
a'p = ap and work only close to the Feshbach resonance. 
Of course this issue deserves further investigations which 
are beyond the scope of the present paper. 



III. HOMOGENEOUS FERMI-BOSE MIXTURE 

We consider a dilute mixture of Np bosons and Np 
fcrmions in a box of volume V. As previosly stated, the 
fermions are equally distributed in two hyperfine states 
and can be considered as a single superfluid Fermi system 
whose energy density is given by Eq. ([1]). Here and in 
the following we will use a local density approximation 
and assume that the interaction parameter y entering the 
universal function f(y) depends on the renormalized s- 
wave scattering length between fermions y = (kpa'p) -1 . 
In our treatment of the BCS-BEC crossover we use a' F 
given by Eq. ([4]) for ap < and a' F = ap for ap > 0. 
The energy density of a uniform, homogeneously mixed 
phase of condensate bosons and superfluid fermions can 
then be written as 
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where (\/2)gBn 2 B is the energy density of a Bose-Einstein 
condensate while gBF^B^F is the contribution from the 
Bose-Fermi interaction, ap and app are not signifi- 
cantly modified by the medium if the dilutcness con- 

1/3 1/3 

ditions asrip <C 1 and aBFn F <C 1 arc satisfied 
[H, HI, Bill, therefore we use for them their bare 
values. 

The mixed phase is energetically stable if its energy is 
a minimum with respect to small variations of the densi- 
ties np and rip, while the total number of fermions and 
bosons are held fixed. To get the equilibrium densities 
one must then minimize the function 



£ = S — fiF np — hb tib 



(7) 



where [If and fj,p are Lagrange multipliers (imposing 
that the numbers of fermions and bosons are fixed) which 
may be identified with the Fermi and Bose chemical po- 
tentials. Setting the derivatives of £ with respect to tif 



and np equal to zero, one finds: 



p F = An 2 J 3 ( f(y) - lf\y) ) + 9bf > (M 



— gBn B + gBF n F 
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The solution of Eqs. © and ©, gives a minimum if the 
corresponding Hessian of £ is positive, i.e. if: 
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implying: 



np < 



Agp 
Sbf 



f{y) 7 yf{y) + y Q y 2 f"(y) 



(ii) 

(Notice that n p appears also in the right side of the in- 
equality via y = ({2>ir 2 np) l l 3 a'p) -1 ) . The solution of this 
inequality gives the region in the parameters' space where 
the homogeneous mixed phase is dynamically stable. Dy- 
namical stability, corresponding to a local minimum of £ 
is a necessary but not sufficient condition for the ener- 
getic stability (also called thermodynamical stability at 
non-zero temperatures) which requires instead the global 
minimum of £ . In fact, we shall show that there could 
exist inhomogeneous, two-phase, configurations with en- 
ergy lower than that of the uniform, one-phase, config- 
uration considered in Eq. ([5]). In Fig. 2 we plot the 
region of dynamical stability of the homogeneous mix- 
ture in the plane (l/a' F ,np), where a' F is given by Eq. 
for ap < and a' F = aj? for ap > 0. The solid line 
corresponds to the parametrization of Eq. © for f(y) in 
solving the equality in pip , while the dashed line follows 
fromEq. ©. 

Fig. 2 shows that in the BCS regime (l/a' F <C —1) the 
critical density np is close to (2/3) 3 . Near the unitarity 
limit (1/cl'f = 1/clf = 0) the critical density strongly re- 
duces and the homogeneous mixture is no longer stable 
for sufficiently large l/a' F > 0. This is not surprising 
because for a positive a' F , where a' F = ap, the Fermi 
component becomes a gas of Bose-condensed molecules 
of mass 2niF, which likes to separate from the atoms of 
mass Til p • As shown in Fig. 2, there exists an upper 
critical value of l/a' F for the stability of the homoge- 
neous mixture. The existence of such upper critical value 
can be explained by the following simple analytical argu- 
ment: in the deep BEC regime one treats a homogeneous 
Bose-Bose mixture as composed of atoms with mass tub 
and molecules with mass uim = Imp and interaction 
strength <?m = knTi 2 aM /mM with au — 0-6a' F but a' F — 
ap. Then, from the stability condition g BF < <7b5m/4, 
one finds the critical value at vanishingly small fcrmion 
density as l/a F = l/a! F = 7r/(3Tr 2 ) 2 ^(Ag B /g 2 B F ) = 
0.20(A<7b/<7bf)- This is exactly the value where both 
solid and dashed lines of Fig. 2 meet the np = axis. 
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FIG. 2: Region of dynamical stability of the homogeneous 
mixture in the plane (l/a^r, 71f). Solid line: obtained with 
the universal function f(y) derived by Kim and Zubarev [l5j . 
Dashed line: obtained with the universal function f(y) de- 
rived by Manini and Salasnich [l6j]. In the inset there is 
a zoom of the region with positive Fermi-Fermi scattering 
length a' F , where a' F — a F - The fermionic density n F is in 
units of A 3 g%/ g% F and the scattering length a' F is in units of 



the homogeneous mixture becomes dynamically unstable 
when the lower branch becomes imaginary. 



IV. TWO-PHASE FERMI-BOSE MIXTURE 

The Bose and Fermi components can form distinct 
phases, which we label by the index i. If we ignore in- 
terpenetration effects, a possible phase-separated config- 
uration is described by the number / of phases present, 
the bosonic densities ng ,, the fermionic densities npi in 
each phase, and the fractions vi of the total volume they 
occupy. Since the total number of particles is given, the 
following relations must hold: 
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(14) 



and Y^i v i = 1- When J = 1 we have the case of a 
homogeneous mixture discussed in the previous section. 
The other possiblity for the system we are considering is 
1 — 2 [47]. In this case the total energy is the sum of 
the contributions due to fermions, to bosons and to their 
mutual interaction: 



Moreover, the asymptotic behavior of the universal func- 
tion f(y) ~ 1/y + 0(l/y 5 / 3 ) for large y, implies that the 
boundary of the stability region in Fig. 2 reaches the axis 
rif = with a positive slope. As a consequence there 
must be an interval of l/a' F values where the stability 
region has both a lower and an upper bound in np . Nu- 
merically, one finds that this reentrant behaviour of the 
stability line occurs for values between 0.2Q(Ag B / g BF ) 
and 0.28(Ag B /g BF ) by using Eq. ^ (see the dashed 
line in the inset of Fig. 2)) and between 0.2Q(Ag B / g BF ) 
and 0.22(Ag B /g BF ), but with much smaller Fermi den- 
sities using Eq. (see the solid line). 

It is interesting to observe that the inequality (fTTj) can 
be written as 



np < 



9Bm F n2 
9 2 bf 



(12) 



where c F = v F (f(y) - \yf(y) + j-,y 2 f" (y)) 1/2 /y/Z is 
the sound velocity of the superfluid Fermi component, 
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with fi = f(yi) and yi = ((37r 2 n Fjl ) 1 / 3 a^) \ where a' F 
is given for ap < by Eq. ([4]) properly renormalizing 
the scattering length between fermions in regions where 
there are bosons. For of > we set a' F = ap. Note that 
we have omitted contributions from interface effects since 
they should be negligible if the system is in a large box 
[48| . Equilibrium requires the equality of the pressures 
in each phase: 

p i = 1 An F.i (fi - ^Viflj + \dB n% ti + g B F n B ,i n F .i , 

where f[ — f'(yi) , and moreover the equality of the 
chemical potentials of each species: 



fi F ,i = AnJ,^ ( fi - ^Vif\, ) + 9bf n B ,i ■ (17) 



with v F = ^2An 2 / 3 /m F the Fermi velocity. The sound for fermions and: 
velocity cbf of the Fermi-Bose mixture can be easily ob- 
tained following a procedure similar to the one suggested ' 
in Ref. [461 ] for a two-component Bose-Einstein conden- 
sate. One finds: 



9b n B s + gBF riF,i 



(18) 



cbf 
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(13) 

where c B = \/gBn<B/ m B is the sound velocity of the 
Bose gas. Thus the sound velocity has two branches and 



for bosons. If the boson density n B ^ is non-zero in both 
phases, then the chemical potentials [i B ^ must be equal. 
If the density of one species vanishes in one phase, then 
its chemical potential in that phase must be higher than 
in the other. For example, if the boson density is zero 
in one phase, then the boson chemical potential in that 
phase must be higher than in the other one, i.e. 



Hb,i > Ms,3-t if n B , 







(19) 
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The same is of course true for fermions. 

Thus, one has to distinguish four cases of equilibrium, 
which must be analyzed one by one: 

(i) Two pure phases: The bosons and fermions are com- 
pletely separated corresponding to n F \ = 0, ns.2 = 
and Ub.i 0, n F ^ 0. 

(ii) A mixed phase and a purely fermionic one: The boson 
density vanishes in one region corresponding to rip^i ^= 0, 
riB,2 = and Ub,x 0, np,2 ^ 0. 

(iii) A mixed phase and a purely bosonic one: The 
fermion density vanishes in one region corresponding to 
the conditions np,i — 0, ns.2 and ns.i 7^ 0, np,2 ^ 0. 

(iv) Two mixed phases: All densities ns,i, tif,i and ub,2, 
np 2 arc different from zero, while ns.i ^ ns.2 and uf.i ^ 
np^2- 

The analysis of these four cases is similar to that 
performed for the Fermi-Bose mixture with a single 
fermionic hyperfine state by Viverit, Pethick and Smith 
[28[ and also by Das 30], who investigated the one- 
dimensional mixture. In our problem the presence of the 
universal function f(y) increases the complexity of the 
algebric manipulations. 

We leave out the algebra and discuss the results as a 
function of the inverse scattering length. First of all, the 
analysis shows that the cases (iii) and (iv) are not realiz- 
able. Thus, a separation into two phases, each with dif- 
ferent, non-zero, concentrations of bosons and fermions 
is never in equilibrium, nor is one with a purely bosonic 
phase and a mixed one. The cases (i) and (ii) are instead 
possible if g B F > 0, while for gsF < these configura- 
tions are unstable. 

In Fig. 3 we plot the phase diagrams (nB,np) of 
the Fermi-Bose mixture for three values of the bare s- 
wave scattering length ap and qbf > 0. In the fig- 
ure there are the results obtained by using Eq. ((3|) to 
model the universal function f(y). Very similar results 
are obtained by using Eq. (2]| . Above the dot-dashed line 
there are two pure phases where bosons and fermions are 
completely demixed and therefore a' F = ap. Between 
the dot-dashed line and the solid line there is a mixed 
phase with a fermion scattering length renormalized by 
the presence of the bosons and a purely fermionic one 
where a' F = ap. Below the solid line bosons and fermions 
are homogeneously mixed. Fig. 3 shows that, as X/ap 
is increased from a negative value to the Feshbach reso- 
nance (1/clf = X/a'p = 0), the region of pure separation 
increases, while the other two are strongly reduced. This 
effect is stronger above the Feshbach resonance (1/ap > 
and a' F = ap), where pairs of Fermi atoms become Bose- 
condensed molecules. As expected from the dynamical 
stability analysis of the previous section, when 1/ap ex- 
ceedes a critical value the region of pure separation oc- 
cupies the whole phase diagram. This critical value is 
1/ap = Q.2%{Ags/ 9 B f) according to Eq. © and it is 
instead 1/ap = 0.22(Ag B / g BF ) using the Eq. @. 

As discussed in the previous section, a physical conse- 
quence of the asymptotic behavior of f(y) for y — > +oo 
is the existence of an interval of 1/of, where the re- 
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FIG. 3: (color online). Phase diagrams in the plane (tib,^f) 
of the Fermi-Bose mixture with gsF > and three values of 
the inverse Fermi- Fermi scattering length «f. Note the differ- 
ent scales of the top panel. Above the dot-dashed line: bosons 
and fermions are completely separated (light grey region). Be- 
tween the dot-dashed line and the solid line: a mixed phase 
and a purely fermionic one. Below the solid line: bosons and 
fermions are fully mixed (dark grey region). The dotted line 
gives the curve below which the Eq. can be used (mixture 
of 6 Li and 23 Na atoms). For a F > we set a' F — a F . The 
bosonic density ub is in units of A Li g%/\gBF\ 5 ■ The fermionic 
density n F is in units of A 3 g%/g% F . The Fermi- Fermi scat- 
tering length a F is in units of g% F / (Ag F )- 
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FIG. 4: (color online). Phase diagram with gsF > and 
1/ap =0.25. Above the dot-dashed line: bosons and fermions 
are completely separated (light grey region). Between the dot- 
dashed line and the solid line: a mixed phase and a purely 
fermionic one coexist. Only in the dark grey region bosons 
and fermions are homogeneously mixed. In the hatched re- 
gions none of the possible configurations is dynamically and 
energetically stable. Units as in Fig. 3. 
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gion of dynamical stability of the mixed phase has not 
only an upper bound in the fermion density but also 
a lower bound. By using Eq. ([3|) this happens for 
0.20 (Ag B /g% F ) < l/a F < 0.28 (Ag B /g 2 BF ). In Fig. 4 
we plot the phase diagram obtained with Eq. ([3]) and 
1/a.F = 0.25. As in Fig. 3, the dark grey region of 
the diagram is where the fully mixed configuration is dy- 
namically and energetically stable. In Fig. 4 there are 
domains, indicated as hatched regions, where none of the 
possible configurations is dynamically and energetically 
stable. Obviously, in these domains the total energy is 
no longer given by Eq. (fTS"]) and the role of the interface 
may be relevant. There is a similar effect for qbf < 0. 
In fact, for gsF < only the single homogeneous phase, 
where bosons and fermions are fully mixed, is stable: in 
Fig. 3 and in Fig. 4, this region is in between the two 
dashed lines, and it is independent on ng. In the other 
regions of the phase diagram all the configurations are 
unstable, and presumibly this corresponds to a collapse 
where fermions and bosons are clumped together [28| . 

Our predictions can be verified by using alkali-metal 
atoms. For the Fermi-Bose mixture of 6 Li atoms and 
23 Na atoms one has clbf = 30 ao, where ao = 0.53xl0~ 10 
m is the Bohr radius [43| . In this case it is easy to find 
that the scaling units of nj? and tib, i.e. (Ags / Qbf) 3 
and A 3 g B /\g B F\ 5 , are both about 10 s fj,m~ 3 . For the 
Fermi-Bose mixture of 40 K atoms and Rb atoms one 
has instead o,bf — — 284 ao [35| and the scaling units are 
about 10 2 /Ltm -3 for fermions and about 10 3 /zm -3 for 
bosons. From these densities and Fig. 3 it follows that 
in the unitarity limit (1/a.F = 0) the mixing-demixing 
transition can be obtained for instance in a mixture of 
~ 10 6 6 Li atoms and ~ 10 6 23 Na atoms in a volume of 
~ 0.5 /im 3 . These numbers of atoms are already realized 
in experiments with ultra-cold alkali- metal atoms. 



V. INCLUSION OF AN EXTERNAL 
POTENTIAL 

In actual experiments, the ultra-cold atoms are usually 
trapped by an external harmonic potential. For simplic- 
ity we consider the same spherical harmonic potential 



U(r) = \k 



(20) 



acting on fermions and bosons, where K = m F ui F = 
rriB^B ^ S t ne elastic constant and uj f and lub are the cor- 
responding harmonic frequencies. The energy functional 
of the mixture can be written as 



E 



-£ 



(n F (r), n B {r)) + U{r) n F (r) + U(r) n B (r)\d 3 r (21) 



where £ (n F (r), ub (r)) is given by Eq. © while the gra- 
dient terms give the quantum pressure for bosons and 



the von Weiszacker contribution to the kinetic energy of 
fermions with a nonuniform density. By minimizing the 
energy with a fixed number Nb of bosons and Nf of 
fermions we get 



flp 



n 2 v 2 y^7(0 

2m F ^fn F (r) 



+AriF(rp 3 (j(y(r)) - ^f{y(r))j +g B Fn B (r)+U(r), 

(22) 



1-lB 



h 2 \7 2 Jn B {r) . , . . TT . . 

7, 1 ? \ +9BnB(r)+gBFnF{r) + U(r) , 

2m B ^n B (r) 

(23) 

where y(r) — ((37r 2 riF(r)) 1 / 3 a F ) x . Again a' F = of in 
the case of full demixing and also for of > 0. The pres- 
ence of a harmonic external potential U(r) strongly mod- 
ifies the conditions for demixing. Clearly, due to the har- 
monic confinement partial demixing is always possible. 
For instance, with gsF = both Fermi and Bose clouds 
occupy the center of the trap but the effective radii can 
be quite different and in this case demixing will appear 
far from the center of the trap. Such conditions can be 
estimated analitically from the previous equations (|22"|) 
and (|23|) within the Thomas-Fermi (TF) approximation, 
i.e. neglecting von Weiszacker and quantum-pressure 
terms. This is a good ap prox imation when gs > 
and Nf and Nb are large [49l [50j . For qbf — the 

TF approximation gives Rb = (ldgsNs / (4Trm,BW B )) 
for the effective radius of the Bose cloud. The ef- 
fective radius Rf of the Fermi cloud depends instead 
on l/a' F . For l/a' F <C —1 (BCS regime) one finds 
R F = (48N F ) 1/6 y/h/(mFUJF); for l/a' F = (unitarity 
limit) one finds R F = /(0) 1 / 4 (48A F ) 1/ 6 ^h/ (m F 0JF); 
for l/a' F = 1/ap ^ 1 (BEC regime) one finds 

Rf = (0.60) 1 / 5 (l5g' F N B /(32Trmpuj 2 F )) 1/5 , where g' F = 
4:Ttha' F /m F . 

Also for gsF > partial demixing is possible by vary- 
ing the number of atoms. This is the typical situation 
of the 6 Li- 23 Na mixture where gsF I 9b = 0.84. In Fig. 
5(a) we plot the density profiles of the Bose and Fermi 
clouds obtained in the unitarity limit l/a' F — with 
N B = 10 6 and N F — 10 5 . In this case the radii of the 
two clouds are very similar but, as shown in Fig. 5(b) 
where Nb — 10 6 and N F = 5 x 10 6 , by increasing the 
number N F of fermions one produces partial demixing 
with Bose and Fermi atoms in the center of the trap and 
an external shell of fermions. The results of Fig. 5 have 
been obtained by inserting the Thomas-Fermi density of 
fermions, Uf{t) 1 into Eq. (|23|) . In this way one gets the 
following nonlinear Schrodinger equation 



ft 2 V 2 
2m b 



+ U(r) + gBriB(r) + gBFnp(r) 



ip(r) 



H B ip(r) 
(24) 
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where tib{t) = IVK 7 *)! 2 an( i VK 7 ") is the macroscopic wave 
function of the Bose condensate. We have solved this 
equation by considering its time-dependent version and 
using a finite-difference Crank-Nicholson method with 
imaginary time (5lj . 

Usually both Fermi and Bose atoms occupy the cen- 
ter of the trap, but under some conditions one of the 
two species can be expelled from the center. In partic- 
ular, as suggested by Molmer [24j], if the bosonic cloud 
is practically independent on the fermions, one may get 
the density profile of bosons as 



(25) 



n B (r) ~ — (jib ~ U(r)) 

9B 



by neglecting the quantum-pressure term in Eq. (|23|) . 
Inserting this result into Eq. (|22p one gets the following 
TF formula for the fermionic density profile nj?(r): 



An F {r)^ (f(y(r)) - ~/'(y(r)) 



9b f TTl s. 9bf 
[is - U[r) 1 

9B V 9B 



(26) 



where the terms proportional to g BF are absent in re- 
gions with vanishing %(r). This formula shows that if 
9bf/9b < 1 the fermions will occupy the center of the 
trap together with bosons. Instead, if g BF 7 1 9b > 1 the 
fermions are expelled from the center of the trap and 
form a shell outside the Bose condensate, i.e there is a 
core of bosons surrounded by a shell of fermions. In Fig. 
5(c) we show the profiles of the 6 Li- 87 Rb mixture with 
the artificial value 9bf I 9b = 5 and setting Nb = 10 6 
and Np = 10 5 . It this case the fermionic cloud is ex- 
pelled from the center of the trap. This effect is more 
clearly shown in Fig. 5(d) where we set N B — 10 6 and 
Np = 5 x 10 6 . Note that by using a much larger value of 
the ratio gBF/gB we find full demixing. 

For the 87 Rb- 40 K mixture the Bose-Fermi scattering 
length is attractive, i.e. Qbf < 0. In this case the Fermi 
atoms are mixed with the Bose atoms up to the collapse. 
As verified by various authors [HI, 13 [El, [53| with a 
Fermi-Bose mixture and with fermions in a unique hyper- 
fine state, the collapse point is well described by the local 
density approximation at the point where the dynamical 
instability sets up. In our case this point is obtained from 
Eq. (HH) as 



n F (0) 



Ag B 
9 2 bf 



f(y c (0))~-y c (0)f'(y c (0)) 



+ I ^ c (0) 2 .f'V(0)) 



(27) 



where y c (0) = ((37r 2 nJ,(0)) 1/3 a' F )- 1 and np(0) is the 
critical density at the center of the trap above which 
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FIG. 5: Fermi-Bose mixture of 6 Li- 23 Na atoms in a harmonic 
trap with a resonant Fermi- Fermi scattering length (1/cif = 
l/a' F = 0). We set as/an = 0.01, where an = '(msws) 
is the characteristic harmonic length of the trap and Nb = 
10 6 . a) N F = 10 5 and g B F / 1 Qb = 0.84; b) N F = 5 x 10 6 and 
g BF /gB = 0.84; c) N F = 10 5 and g BF /g B = 5; d) N F = 
5 x 10 6 and gBF/gB = 5. Solid line: density profile np(r) of 
the cloud of two-spin 6 Li atoms; dashed line: density profile 
n B (r) of the cloud of 23 Na atoms. Lengths are in units of an 
and densities in units of aT? . 



there is the collapse. The plot of n c F (0) as a function 
of l/a' F is precisely the solid (dashed) curve of Fig. 2 
taking f(y) from Eq. © (Eq. ©). Using the most 
recent available data on the 40 K- 87 Rb mixture [35| we 
predict that in the deep BCS regime (l/a' F -C — 1) the 
critical fermionic density is n c F (0) — 140 //m -3 , while in 
the unitarity limit (l/a' F = 0) it reduces to n F (0) = 10 
/im~ 3 . 



VI. CONCLUSIONS 

We have investigated a Fermi-Bose mixture with the 
Fermi atoms in two equally populated hyperfine states. 
In particular, we have analyzed the mixture as the Fermi- 
Fermi scattering length is varied across a Feshbach reso- 
nance by using two efficient parametrizations of the uni- 
versal function f(y), which characterizes the BCS-BEC 
crossover. We have found that the phase diagram of the 
system strongly depends on the Fermi-Fermi scattering 
length. The conditions for demixing and collapse have 
been studied in a box and also in a harmonic trap by 
using a density functional approach. Our results sug- 
gest that in the unitarity limit, where the renormalized 
scattering length of the Fermi-Fermi interaction goes to 
infinity, demixing or partial demixing can be observed 
using a mixture of 6 Li and 23 Na atoms. Collapse can be 
instead obtained using a mixture of 40 K and 87 Rb atoms 
and varying the fermionic density of the system. Many 
interesting issues regarding Fermi-Bose mixtures across a 
Feshbach resonance remain to be investigated both the- 
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